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FOREWORD 


This  document  is  a  final  task  report  for  the  RSVP  Upgrade  md  Mainten¬ 
ance  Task  completed  under  Contract  No.  N00014-7&-C-0818  for  the  Naval  Ocean 
Research  and  Development  Activity. 

The  RSVP  system  retrieves  sound  speed  profiles  from  a  modified  NODC 
(National  Oceanographic  Data  Center)  data  base,  and  determines  representative 
profiles  from  the  retrieved  profiles.  These  representative  profiles  are  then  used 
for  oceanographic  analysis,  including  acoustic  modeling. 

Under  this  task,  modifications  were  made  to  the  computer  programs  in  the 
RSVP  system  to  make  them  easier  to  use.  This  task  also  included  conversion  to 
new  computer  facilities,  data  base  maintenance,  and  documentation  of  an  exist¬ 
ing  undocumented  program. 

This  document  is  in  three  parts.  The  first  part  describes  the  improvements 
made  to  the  RSVP  retrieval  program.  The  second  part  is  the  documentation  of 
the  RSVP  plotting  program.  The  third  part  is  a  description  of  the  Cyber  Control 
Language  procedure  that  is  used  to  run  the  RSVP  system. 

ODSI  thanks  Eigoro  Hashimoto  formerly  of  NORDA  Code  320  and  George 
Kerr  of  NORDA  Code  321  for  assistance  with  program  specifications  and  for 
furnishing  an  environment  conducive  to  completion  of  this  task. 


OCEAN  DATA  SYSTEMS,  INC. 


tk  Dr.  A.  Anderson  (Head  -  NORDA  Code  320) 

rNM:  Jacqueline  Webster 

dm:  30  June  19S0 

lusiMt:  Upgrade  of  the  RSVP  System 

The  RSVP  system  retrieves  sound  profiles  from  a  modified  and  reformatted 
NODC  ocean  station  data  base,  and  determines  representative  profiles  from 
among  those  that  are  retrieved.  The  representatives  are  then  used  as  input  to 
acoustic  modelling  programs  and  for  other  oceanographic  analysis.  Under 
Contract  No.  N00014-78-C-0818,  modifications  were  made  to  two  of  the 
programs  in  the  system,  RSVP  and  RSVPPLP.  These  modifications  are  summar¬ 
ized  as  follows: 

Changes  to  RSVP 

1.  Automation  of  the  procedure  whereby  deep  ocean  sound  speeds  are 
retrieved.  Deep  ocean  sound  speeds  are  now  automatically  retrieved 
from  a  disk  file  rather  than  manual  user  selection  from  a  punched  card 
deck  being  required. 

2.  Extension  of  the  summary  of  soiard  speed  profiles  by  year  from  1974  to 
1999. 

3.  Addition  of  parity  error  checking  on  input  of  sound  speed  profiles. 

4.  Clarification  of  printer  output. 

3.  Addition  of  a  new  format  for  disk  output  (used  as  input  to  program 
RSVPPLP). 


Changes  to  RSVPPLP 

1.  Increasing  the  flexibility  of  input  specifications. 

2.  Addition  of  new  features  to  the  plotter  output  and  making  it  more 
general. 

3.  Uncovering  and  correction  of  program  errors. 

4.  Modification  of  the  card  output  to  standard  FACT  format. 

Attached  is  documentation  of  the  changes  made  to  RSVP  along  with 
Instructions  for  its  use. 


Documentation  of  RSVPPLP,  prepared  as  part  of  the  work  performed  under 
the  previously  mentioned  contract,  appears  as  a  separate  task  report.  RSVP  is 
documented  in: 

AESD  Sound-Speed  Profile  Retrieval  System  (RSVP),  AESD  Technical 
Note  TN-74-03,  October  1974,  3ohn  3.  Audet,  3r.  and  Gregory  G. 
Vega. 


cc:  Eigoro  Hashimoto  -  NOROA  Code  320 
Gilbert  3acobs 


I.  Changes  to  program  RSVP 
1.1  Archival  deep  sound  speeds 

The  RSVP  system  uses  archival  deep  sound  speeds  to  determine  critical 
depth  where  available  data  are  not  sufficiently  deep,  and  to  extend  sound  speed 
profiles  when  specified  by  the  user.  Previously,  the  RSVP  user  had  to  find  the 
appropriate  card  containing  the  deep  sound  speeds  from  a  card  file,  duplicate  it 
and  manually  insert  it  in  his  input  deck.  The  program  has  been  modifi^  so  that 
the  deep  sound  speeds  can  be  specified  in  one  of  three  ways: 

1.  The  deep  sound  speeds  can  be  specified  on  a  punched  card,  as  before. 

2.  The  Marsden  square  number  and  quadrant  are  calculated  by  the  pro¬ 
gram,  and  then  used  to  search  for  matching  deep  sound  speeds  in  a  disk 
file. 

3.  The  Marsden  square  number  and  quadrant  can  be  specified  by  the  user 
on  punched  card  input.  These  are  then  used  to  obtain  the  matching 
deep  sound  speeds  from  the  disk  file. 

The  file  containing  the  archival  deep  sound  speeds  consists  of  one  record 
for  each  Marsden  square  quadrant.  A  description  of  the  record  follows: 


Format 

Character 

Positions 

Variable 

Name 

Description 

F6.I 

1-6* 

BLT(l) 

Maximum  depth  in  meters  for  the 
I-th  request  area. 

10F6.1 

7-12,* 

13-18, 

19-20, 

• 

(SVEXT(J,I), 

3=1,10) 

Archival  deep  sound  speeds  in  me¬ 
ters/sec  for  the  I-th  request  area. 

• 

6i-66 

13 

69-71** 

LMSQ 

Marsden  square  number. 

11 

73 

LMSQ5 

Quadrant  number. 

*  The  value  is  placed  anywhere  in  the  positions  specified  with  the  decimal  point 
expressed. 

**The  value  is  right-justified  in  the  positions  specified  with  no  decimal  point. 
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1.2  Summary  by  Year 

The  printed  output  from  RSVP  contains  a  count  of  the  number  of  profiles  in 
the  area  of  interest  for  each  5-year  period  between  1900  and  1974.  Added  to 
this  summary  is  a  count  of  the  profiles  between  1975  and  1999.  This  was  done  in 
preparation  for  the  updates  that  will  be  made  to  the  sound  q>eed  data  base. 

1.3  Parity  Errors 


In  the  previous  version  of  RSVP,  a  hardware  error  encountered  on  reading 
the  sound  speeds  was  treated  in  the  same  way  as  an  end-of-data.  Thus  a  bad 
tape  or  disk  could  cause  a  user  to  assume  that  there  was  no  data  in  the  area, 
when  there  actually  was  data  there.  The  current  version  of  RSVP  prints  on  error 
message  and  aborts  the  job. 

1.4  Printer  Output 

The  headings  in  the  printer  output  from  RSVP  were  modified  in  order  to 
make  it  easier  to  see  what  data  were  available  in  the  data  base,  and  what  data 
were  selected  for  analysis.  A  sample  of  the  RSVP  printer  output  is  presented  in 
Section  II.4. 

1.5  Disk  Output 

The  previous  version  of  RSVP  wrote  all  the  printer  output  to  a  disk  file  for 
use  by  the  next  program  in  the  system,  RSVPPLP.  This  file  used  more  disk  space 
than  necessary,  causing  a  shortage  of  scarce  disk  storage  when  the  RSVP  system 
was  used  in  preparing  for  an  exercise.  Therefore,  the  disk  file  output  from  RSVP 
was  changed  to  a  simpler  more  compact  version. 

The  disk  file  consists  of  sets  of  sound  speed  profiles  retrieved  by  RSVP  as 
well  as  information  used  in  retrieval  of  the  profiles.  Each  set  in  this  file  is  the 
output  from  one  request  (by  ocean  location  and  time  of  year)  to  program  RSVP. 
There  is  a  maximum  of  16  sets  per  file. 

Each  set  in  the  file  consists  of  one  record  each  of  Types  1-4,  followed  by 
35  records  of  Type  5  as  described  on  the  following  pages.  Record  Types  1-4  are 
images  of  the  original  request  specification  input  to  program  RSVP.  Records  of 
Type  5  contain  the  sound  speed  profiles  that  were  retrieved  by  RSVP  as  a  result 
of  the  request  shown  in  Records  1-4. 

A  description  of  the  five  record  types  follows. 
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Record 

Character 

Variable 

Type 

Format 

Positions 

Name 

Description 

1 

6  A 10 

1^0 

(REQT(K,l), 

Up  Id  60  alphanumeric  char- 

K*l,6> 

acters  for  an  Identifying  title 
and  a  date  lor  the  1-th  request. 

12 

78-79** 

I 

Request  number  (1  -  16). 

11 

80 

ISEQ 

Record  type  number*  equal  to 

1. 

2 

12 

1-2** 

Ltd, I) 

Degrees  of  lower  latitute 

12 

4-5** 

LL<2,I) 

Minutes  limit  of  search 

A1 

7 

LL(3,1) 

Hemisphere  area  lor  1-th  re- 

quest. 

12 

10-11** 

LL<4,1) 

Degrees  of  upper  latitude 

12 

13-14** 

LL(5,l) 

Minutes  limit  oi  search 

A1 

16 

LL<6,I) 

Hemisphere  area  for  I-th  re- 

quest 

13 

21-23** 

LLC7,1) 

Degrees  ol  lower  longi- 

12 

25-26** 

LL(8,1) 

Minutes  tude  limit  ol 

A1 

28 

LL(9,1> 

Hemisphere  search  area  for 

I-th  request 

13 

31-33** 

LL(10,1) 

Degrees  of  upper  longi- 

12 

35-36** 

LL(11,1) 

Minutes  tude  limit  ol 

1 

1 

A1 

38 

LL(12,1) 

Hemisphere  search  area  for 

! 

I-th  request. 

■ 

13 

50-52** 

LMSQ 

Marsden  square  number  and 

I 

1 

11 

54 

LQUAD 

quadrant  number  used  by  RSVP 
to  look  up  the  deep  sound 
^eds  in  the  archival  file.  11 
blank*  deep  sound  speeds  are 
specified  by  INDM  below. 

A1 

55 

INDM 

If  LMSQ  and  LQUAD  above  are 
specified*  INDM  is  disregarded. 

U  INDM  is  ”M",  the  Marsden 
square  and  quadrant  are  cal- 

culated  from  latitude  and  longi¬ 
tude.  If  INDM  is  blank,  deep 
sound  speeds  are  entered  from 
cards  of  type  4. 

12 

78-79** 

1 

Request  number  (1-16) 

11 

80 

ISEQ 

Record  type  number,  equal  to 

2. 

*The  vaiu«  it  placed  anywhere  in  the  positiont  4>ecified  with  the  decimal  point 
expretted. 

**The  value  it  right>)uttliied  in  the  positiont  tpecUied  with  no  decimal  point. 
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Record 

Type 

Format 

Character 

Positions 

Variable 

Name 

Description 

3 

12(I2,1X) 

1-2,** 

4-5, 

7-8, 

: 

(1SN(M,I), 

Msi,12) 

Month  numbers  (01-12)  used  in 
data  search  for  1-th  request. 

34‘-35 

14 

38-41** 

IY(l) 

Earliest  year  used  in  data 
search  for  1-th  request. 

14 

44-47** 

MIND(I) 

Minimum  search  depth  in  me¬ 
ters  for  I-th  request. 

F4.0 

50-53* 

XSONiCd) 

Maximum  search  depth  in  me¬ 
ters  of  sonic  layer  for  I-th  re¬ 
quest. 

F4.0 

56-59* 

SSXMlNd) 

Minimum  search  depth  in  me¬ 
ters  of  subsurface  maximum  for 
1-th  request. 

12 

78-79** 

1 

Request  number  (1-16) 

11 

80 

ISEQ 

Record  type  number,  equal  to 

3. 

4 

11F6,1 

1-6* 

BLTd) 

Maximum  bottom  depth  in  me¬ 
ters  for  the  Marsden  square 
quadrant. 

7-i2,* 

13-18, 

m 

• 

61-66 

(SVEXT(3,I), 

3=1,10) 

Archival  deep  sound  speeds  in 
meters/sec.  used  to  extend  pro¬ 
files  for  the  1-th  request. 

69-71** 

MSQ 

Marsden  square  number  lor  ar¬ 
chival  deep  sound  speeds. 

n 

73 

MSQ5 

Marsden  square  quandrant. 

12 

78-79** 

I 

Request  number  (1-16) 

n 

80 

ISEQ 

Record  type  number,  equal  to 

4. 

There  is  one  record  oi  each  oi  Types  representing  an  RSVP  reqviest. 

*The  value  is  placed  anywhere  In  the  posltiorts  specif  led  with  the  decimal  point 
expressed. 

**The  value  is  right>justified  in  the  positions  specified  with  no  decimal  point. 
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Record 

Type 

Format 

Character 

Positions 

Variable 

Name 

Description 

5 

12 

9-10« 

3 

Sequence  number  (1-35)  for 
sound  speed  profile  point. 

15 

16-20** 

LDEPTHO) 

Depth  in  meters  at  3-th  sound 
speed  profile  point. 

F6.1 

25-30* 

RPROFD) 

Sound  speed  in  meters/sec.  at 
the  3-th  depth  for  the  model 
(most  typical)  profile. 

F6.1 

35-40* 

SVLMO,!) 

Mean  sound  speed  in  me¬ 
ters/sec.  at  the  3-th  depth,  for 
all  profiles  selected  for  the  I-th 
request. 

F6.1 

44-49* 

SVMIN(3,I) 

Minimum  sound  speed  in  me¬ 
ters/sec.  at  the  3-th  depth,  for 
all  profiles  selected  for  the  1-th 
request. 

F6.1 

53-58* 

SVMAX(J,I) 

Maximum  sound  speed  in  me¬ 
ters/sec.  at  the  3-th  depth,  for 
all  profiles  selected  for  the  1-th 
request. 

F4.1 

67-70* 

SVLSD(3,I) 

Standard  deviation  from  the 
mean  in  meters/sec.  at  the  3-th 
depth,  for  all  profiles  selected 
for  the  I-th  request. 

15 

81-85** 

CNT(3,I) 

Number  of  observations  at  the 
3-th  depth  for  all  profiles 
selected  for  the  1-th  request. 

F6.1 

95-100,* 

(SVL3(3,K), 

Sotmd  speed  in  meters/sec.  at 

F6.1 

F6.1 

104-109,* 

113-118* 

K=l,3) 

the  3-th  depth  for  the  3  most 
representative  profiles  selected 
for  the  I-th  request. 

15 

123-127** 

LDEPTH(3) 

Depth  in  meters  at  3-th  sound 
speed  profile  point  (repeat  of 
second  field  above) 

*The  value  is  placed  anywhere  in  the  positions  specified  with  the  decimal  point 
expressed. 


««The  value  is  right- justified  in  the  positions  specified  with  no  decimal  point. 
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n.  RSVP  User's  Guide 
n.l  File  Identification 
n.l.l.  Program  Files 

The  FORTRAN  source  oi  program  RSVP  is  stored  in  the  Program  Library 
file  PVEHRSVPPL.  The  binary  form  of  RSVP  is  stored  as  file  PVEHRSVPLGO. 
The  ID  at  Eglin  is  KM,  and  the  ID  at  NSRDC  is  PVEH. 

II.  1.2  Data  Files 

The  archival  deep  sound  speeds  are  stored  on  a  card-image  disk  file  with 
the  name  PVEHDEEPSSARCHIVAL.  The  ID  at  Eglin  is  KM  and  the  ID  at  NSRDC 
is  PVEH. 

At  Eglin  the  sound  speed  data  base  is  stored  on  magnetic  tape  (9-track, 
labelled,  density  SOO  bpi).  At  NSRDC  the  data  base  is  stor^  on  removable  disk 
packs.  The  data  are  grouped  by  geographical  area.  Following  is  a  table  showing 
file  names,  tape  reel  identifiers,  and  disk  pack  identifiers  for  each  group  of  data. 


STORAGE  MEDIA  FOR  THE  RSVP  DATA  BASE 


The  external  identifiers  listed  here  are  not  the  same  as  the  internal  labels  (VSN)  because  ttese  tapes  were 
created  at  NSRDC.  Therefore  the  LABEL  card  used  to  specify  the  tape  should  not  contain  a  VSN. 


created  at  NSRDC.  Therefore  the  LABEL  card  used  to  specify  the  tape  should  not  contain  a  VSN. 


The  external  identifiers  listed  here  are  not  the  same  as  the  internal  labels  (VSN)  because  these  tapes 
created  at  NSRDC.  Therefore  the  LABEL  card  used  to  specify  the  tape  should  not  contain  a  VSN. 


n.2  RSVP  Card  Input 


A  maximum  of  16  ocean  areas  in  whidi  it  is  desired  to  find  typical  somd 
speed  profiles  can  be  ^ecified  in  one  RSVP  run.  For  each  ocean  area 
specification,  referred  to  as  a  request,  there  are  up  to  4  input  cards.  Dt4>licate 
input  cards  may  be  omitted  when  more  than  one  request  is  being  specified.  That 
is,  if  an  input  card  in  one  request  is  the  same  as  the  corresponding  type  of  card 
for  the  previous  request  in  the  run,  then  that  card  may  be  omitted.  A 
description  of  the  4  types  of  input  cards  to  RSVP  follows.  Sample  RSVP  input 
cards  appear  following  the  card  input  description. 
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Card 

Type 

Format 

Character 

Positions 

Variable 

Name 

Description 

1 

6  A 10 

1-60 

(REQT{K,a 

K=l,6) 

Up  to  60  alphanumeric  char¬ 
acters  for  an  identifying  title 
and  a  date  for  the  1-th  request. 

12 

78-79** 

1 

Request  number  (1  -  16). 

11 

80 

ISEQ 

Card  type  numberi  equal  to  1. 

2 

12 

1-2** 

LL(1,1) 

Degrees  of  lower  latitute 

12 

4-5** 

LL(24) 

Minutes  limit  of  search 

A1 

7 

LL(3^) 

Hemisphere  area  for  1-th  re¬ 
quest. 

12 

10-11** 

LL(4,I) 

Degrees  of  upper  latitude 

12 

13-14** 

LL(5,1) 

Minutes  limit  of  search 

AI 

16 

LL(6,I) 

Hemisphere  area  for  1-th  re¬ 
quest 

13 

21-23** 

LL(7.I) 

Degrees  of  lower  longi- 

12 

25-26** 

LL(8,1) 

Minutes  tude  limit  of 

Al 

28 

LL(9,I) 

Hemisphere  search  area  for 

1-th  request 

13 

31-33** 

LL(10,I) 

Degrees  of  upper  longi- 

12 

35-36** 

LL(11,I) 

Minutes  tude  limit  of 

Al 

38 

LL(12,1) 

Hemisphere  search  area  for 

I-th  request. 

13 

50-52** 

LMSQ 

Marsden  square  number  and 

11 

54 

LQUAD 

quadrant  number  used  by  RSVP 
to  look  ip  the  deep  sound 
speeds  in  the  archival  file.  If 
blank,  deep  sound  speeds  are 
^ecified  by  INDM  below. 

Al 

55 

INDM 

If  LMSQ  and  LQUAD  above  are 
specified,  INDM  is  disregarded. 

If  INDM  is  "M”,  the  Marsden 
square  and  quadrant  are  cal¬ 
culated  from  latitude  and  longi¬ 
tude.  If  INDM  is  blank,  deep 
sound  q>eeds  are  entered  from 
cards  of  type  4. 

12 

78-79** 

1 

Request  ruimber  (1-16) 

11 

80 

ISEQ 

Card  type  number,  equal  to  2. 

*The  value  is  placed  anywhere  in  the  positions  specified  with  the  decimal  point 
expressed 


**The  value  is  right-justified  in  the  positions  specified  with  no  decimal  point. 


Card 

Type 

Format 

Character 

Positions 

Variable 

Name 

Description 

3 

1202,1  X) 

1-2,** 

<►-5, 

7-8, 

• 

• 

34*-35 

QSN(M,I), 

M=l,12) 

Month  numbers  (01-12)  used  In 
data  search  for  I-th  request. 

14 

38-41** 

IY(I) 

Earliest  year  used  in  data 
search  for  I-th  request. 

14 

44-47** 

MINDO) 

Minimum  search  depth  in  me¬ 
ters  for  I-th  request. 

F4.0 

50-53* 

XSONICO) 

Maximum  search  depth  in  me¬ 
ters  of  sonic  layer  for  I-th  re¬ 
quest. 

F4.0 

56-59* 

SSXMINO) 

Minimum  search  depth  in  me¬ 
ters  of  subsurface  maximum  for 
I-th  request. 

12 

78-79** 

I 

Request  number  (1-16) 

11 

80 

ISEQ 

Card  type  number,  equal  to  3. 

Cards  of  type  4  are  omitted  when  the  archival  deep  sound  speeds  are  to  be 
obtained  from  the  disk  file,  that  is,  when  columns  50-55  on  the  type  2  card  are 
blank. 


4 

11F6.1 

1-6* 

BLT(I) 

Maximum  bottom  depth  in  me¬ 
ters  for  the  Marsden  square 
quadrant. 

7-12,* 

13-18, 

• 

m 

(SVEXT(3,I), 

3=1,10) 

Archival  deep  sound  speeds  in 
meters/sec.  used  to  extend  pro¬ 
files  for  the  I-th  request. 

61-66 

13 

69-71** 

MSQ 

Marsden  square  number  for  ar¬ 
chival  deep  sound  speeds. 

11 

73 

MSQ5 

Marsden  square  quandrant. 

12 

78-79** 

I 

Request  number  (1-16) 

11 

80 

ISEQ 

Card  type  number,  equal  to  4. 

*The  value  is  placed  anywhere  In  the  positions  specified  with  the  decimal  point 
expressed. 

**The  value  it  right-justified  in  the  positions  specified  with  no  decimal  point. 
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n.3  Sample  jobs 

n.3.1  Sample  Job  Deck  For  Eglin 


Xl07tT37S*NTltP<i.  HASHlMOtO  1609  921BVT16H021  26 

REOUEST  t T APt60  « •PF « SN*USC7 1 . 

ATTACH«TAPE4«P</CH0E£P6SAPCHlVAt.tl0«KH»SN«USETl«MRal. 
LAPELfTAPEl«LaPUJAHEl)SEA«RtU«HO.  PHYSICAL  REEL  NO.  0572 
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FOREWORD 


This  document  is  a  task  report  prepared  under  Contract  No.  N00014-78-C- 
0818  for  the  Office  of  Naval  Research,  Washington,  D.C.,  in  support  of  the  Naval 
Ocean  Research  and  Development  Activity  (NORDA)  >  Code  320,  Bay  St.  Louis, 
Mississippi.  The  report  describes  the  work  performed  in  fulfilling  part  of  the 
contract,  and  serves  as  documentation  for  a  computer  program  to  be  used  in 
plotting  or  punching  representative  sound  speed  profiles  from  an  archival  file  of 
ocean  station  data. 

Ocean  Data  Systems,  Inc.  is  indebted  to  Eigoro  Hashimoto,  NORDA  Code 
320,  for  assistance  with  program  specifications  and  for  furnishing  an  environ¬ 
ment  conducive  to  completion  of  this  task. 
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I.  INTRODUCTION 


The  RSVP  system  was  developed  by  the  Acoustic  Environmental  Support 
Detachment  (AESD)  of  the  Office  of  Naval  Research.  It  retrieves,  along  with 
other  relevant  data  for  acoustic  modeling,  representative  or  modei  sound  speed 
profiles  for  a  given  ocean  location  and  time  of  year.  The  data  base  used  by  the 
system  is  a  si^set  of  ocean  station  data  from  the  National  Oceanographic  Data 
Center  (NODC).  For  further  information  about  the  retrieval  process,  see  the 
following  report: 

"AESD  Sound-Speed  Profile  Retrieval  System  (RSVP),"  by  3ohn  3.  Audet, 

3r.,  and  Gregory  G.  Vega;  AESD  Technical  Note  TN-7<f-03,  October  1974. 

The  report  herein  presented  describes  RSVPPLP,  a  program  for  plotting, 
modifying,  filtering  or  punching  the  RSVP  retrieved  sound  speed  profiles. 
RSVPPLP  is  a  modification  of  program  RSVPPREP,  which  was  part  of  the  RSVP 
system  developed  by  Vega  and  Audet.  Figure  I-l  on  the  following  page  shows 
hiow  RSVPPLP  fits  into  the  RSVP  system. 

Section  11  of  this  document  describes  the  input  required  by  program 
RSVPPLP.  Section  111  describes  the  output  produced  by  the  program.  Section  IV 
contains  individual  descriptions  of  each  module  of  RSVPPLP,  along  with  a 
description  of  COMMON  storage.  Section  V  has  cross  reference  charts  for 
COMMON  storage  and  for  the  program  modules.  Section  VI  contains  information 
needed  for  running  the  program  -  file  names  and  sample  jobs. 
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n.  INPUT 


n.l  Mass  Storage  Input 

Mass  storage  ir>put  for  RSVPPLP  consists  of  a  file  containing  sets  of  sound 
speed  profiles  that  were  output  from  program  RSVP  along  with  information  used 
in  retrieval  of  the  profiles.  Each  set  in  this  file  is  the  output  from  one  request 
(by  ocean  location  and  time  of  year)  to  program  RSVP.  There  is  a  maximum  of 
16  sets  per  file. 

Each  set  in  the  file  consists  of  one  record  each  of  Types  1>4,  followed  by 
35  records  of  Type  5  as  described  on  the  following  pages. 

Record  Types  1-4  are  images  of  the  original  request  specification  input  to 
program  RSVP.  None  of  the  information  rm  Record  Types  1-4,  except  the  title 
from  Record  Type  1  and  the  archival  deep  sound  speeds  from  Record  Type  4,  is 
used  by  program  RSVPPLP.  Records  of  Type  5  contain  the  sound  speed  profiles 
that  were  retrieved  by  RSVP  as  a  result  of  the  request  shown  in  Records  1-4. 

A  description  of  the  five  record  types  follows. 
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taok 


Record 

Character 

Variable 

Type 

Format 

Positions 

Name 

Description 

1 

6A10 

1-60 

(REQT(K,I)b 

Up  to  60  alphanumeric  char- 

K=l,6) 

acters  for  an  identifying  title 
and  a  date  for  the  I-th  request. 

12 

78-79** 

1 

Request  number  (1  -  16). 

11 

80 

ISEQ 

Record  type  number,  equal  to 

1. 

2 

12 

1-2** 

LL(1,1) 

Degrees  of  lower  latitute 

12 

4-5** 

LL(2,1) 

Minutes  limit  of  search 

A1 

7 

LL(3,I) 

Hemisphere  area  for  1-th  re- 

quest. 

12 

10-11** 

LL(4,1) 

Degrees  of  upper  latitude 

12 

13-14** 

LL(5,1) 

Minutes  limit  of  search 

A1 

16 

LL(6,I) 

Hemisphere  area  for  1-th  re- 

quest 

13 

21-23** 

LL(7,I) 

Degrees  of  lower  longi- 

12 

25-26** 

LL(8,1) 

Minutes  tude  limit  of 

A1 

28 

LL(9,1) 

Hemisphere  search  area  for 

1-th  request 

13 

31-33** 

LL(10,I) 

Degrees  of  upper  longi- 

12 

35-36** 

LL(11,I) 

Minutes  tude  limit  of 

A1 

38 

LL(12,1) 

Hemisphere  search  area  for 

1-th  request. 

13 

50-52** 

LMSQ 

Marsden  square  number  and 

54 

LQUAD 

quadrant  number  used  by  RSVP 
to  look  If)  the  deep  sound 
speeds  in  the  archival  file.  If 

blank,  deep  sound  speeds  are 
specified  by  INDM  below. 

A1 

55 

INDM 

If  LMSQ  and  LQUAD  above  are 
specified,  INDM  is  disregarded. 

If  INDM  is  "M",  the  Marsden 
square  and  quadrant  are  cal- 

culated  from  latitude  and  longi¬ 
tude.  If  INDM  is  biank,  deep 
sound  speeds  are  entered  from 
cards  of  type  4. 

12 

78-79** 

I 

Request  number  (1-16) 

11 

80 

ISEQ 

Record  type  number,  equal  to 

2. 

*The  value  is  placed  anywhere  in  the  positions  specified  with  the  decimal  point 
expressed. 


**The  value  is  right-justified  in  the  positions  specified  with  no  decimal  point. 
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Record 

Type 

Format 

Character 

Positions 

Variable 

Name 

Description 

3 

12(12,  iX) 

1-2,** 

7-8, 

• 

(ISN(M,I), 

M=l,12) 

Month  numbers  (01-12)  used  in 
data  seau'ch  for  I-th  request. 

34-35 

14 

38-41** 

1Y(I) 

Earliest  year  used  in  data 
search  for  I-th  request. 

14 

44-47** 

MIND(I) 

Minimum  search  depth  in 
meters  for  I-th  request. 

F4.0 

50-53* 

XSONIC(I) 

Maximum  search  depth  in 
meters  of  sonic  layer  for  1-th 
request. 

F4.0 

56-59* 

SSXMIN(I) 

Minimum  search  depth  in 
meters  of  subsurface  maximum 
for  I-th  request. 

12 

78-79** 

I 

Request  number  (1-16) 

11 

80 

ISEQ 

Record  type  number,  equal  to 

3. 

4 

1IF6.I 

1-6* 

BLTd) 

Maximum  bottom  depth  in  me¬ 
ters  for  the  Marsden  square 
quadrant. 

7-12,* 

13-18, 

• 

• 

61-66 

(SVEXT(3,I), 

3=1,10) 

Archival  deep  sound  speeds  in 
meters/sec.  used  to  extend  pro¬ 
files  for  the  1-th  request. 

13 

69-71** 

MSQ 

Marsden  square  number  for 
archival  deep  sound  speeds. 

11 

73 

MSQ5 

Marsden  square  quandrant. 

12 

78-79** 

I 

Request  number  (1-16) 

11 

80 

ISEQ 

Record  type  number,  equal  to 

4. 

There  Is  one  record  of  each  of  Types  representing  an  RSVP  request. 

*The  value  is  placed  anywhere  in  the  positions  specified  with  the  decimal  point 
expressed. 

**The  value  is  right>justified  in  the  positions  specified  with  no  decimal  point. 
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Record 

Type 

Format 

Character 

Positions 

Variable 

Name 

Description 

5 

12 

9-10** 

2 

Sequence  number  (1-35)  for 
sound  speed  profile  point. 

15 

16-20** 

LDEPTHO) 

Depth  in  meters  at  3-th  sound 
speed  profile  point. 

F6.1 

25-30* 

RPROF(J) 

Sound  speed  in  meters/sec.  at 
the  3-th  depth  for  the  model 
(most  typical)  profile. 

F6.1 

35-40* 

SVLM(3,I) 

Mean  sound  speed  in 

meters/sec.  at  the  3-th  depth, 
for  all  profiles  selected  for  the 

I-th  request. 

F6.1 

44-49* 

SVMIN(3,I) 

Minimum  sound  speed  in 
meters/sec.  at  the  3-th  depth, 
for  all  profiles  selected  for  the 

I-th  request. 

F6.1 

53-58* 

SVMAX(3,I) 

Maximum  sound  speed  in 
meters/sec.  at  the  3-th  depth, 
for  all  profiles  selected  for  the 

1-th  request. 

Fk.l 

67-70* 

SVLSD(3,I) 

Standard  deviation  from  the 
mean  in  meters/sec.  at  the  3-th 
depth;  for  all  profiles  selected 
for  the  I-th  request. 

15 

81-85** 

CNTO,I) 

Number  of  observations  at  the 
3-th  depth  for  all  profiles 
selected  for  the  I-th  request. 

F6.1 

95-100* 

(SVL3(a,K), 

Sound  speed  in  meters/sec.  at 

F6.1 

F6.i 

104-109,* 

113-118* 

K=l,3) 

the  3-th  depth  for  the  3  most 
representative  profiles  selected 
for  the  I-th  request. 

15 

123-127** 

LDEPTHO) 

Depth  in  meters  at  3-th  sound 
speed  profile  point  (repeat  of 
second  field  above) 

*The  value  is  placed  anywhere  in  the  positions  specified  with  the  decimal  point 
e}q>ressed. 

**The  value  is  right*justified  in  the  positions  specified  with  no  decimal  point. 


There  are  35  records  of  Type  5|  one  for  each  profile  depth  level. 

Sets  of  records  of  Types  1-4  followed  by  35  records  of  Type  5  are  repeated  for 
each  RSVP  request,  up  to  a  maximum  of  16  requests. 
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n.2  Card  Input 


Card  input  to  RSVPPLP  consists  of  one  card  for  each  set  of  sound  speed 
profiles  to  be  plotted  and/or  punched  from  an  RSVP  output  file.  The  card 
specifies: 


1.  The  sequence  number  of  the  RSVP  request  to  be  processed  from  the 
RSVP  output  file.  This  corresponds  to  the  number  punched  in  columris  78-79  in 
the  original  request  specification  input  to  the  RSVP  retrieval. 

2.  The  indices  of  the  profiles  to  be  selected  from  those  produced  by  the 
RSVP  request.  Except  for  the  numbers  S  and  6,  the  indices  refer  to  the  column 
numbers  on  the  printed  output  of  the  sound  speed  profiles  from  program  RSVP. 
The  index  numbers  and  the  corresponding  profile  selected  are  shown  below. 

Profile  Index  No.  Profile  Selected 


1 

2 

3 

4 

5 

6 

7 

8 
9 


Model 

Mean 

Minimum 

Maximum 

Mean  minus  standard  deviation 
Mean  plus  standard  deviation 
Representative  No.  1 
Representative  No.  2 
Representative  No.  3 


3.  The  processing  options.  There  are  three  of  these: 

Code  Meaning 

T  Extend  the  specified  profile  with  mean 

values  first  and  then,  if  necessary,  with 
archival  deep  sound  speeds  (from  Record 
Type  4  described  in  Section  II.  1)  to  a 
specified  level. 

M  Extend  with  mean  values  only. 

F  Filter  the  data  by  eliminating  points  from 

the  profile,  but  still  leaving  a  good  repre¬ 

sentation  of  the  curve.  See  the  description 
of  subroutine  FRITZ  in  Section  IV.5  for 
more  details  on  the  filtering  option. 

Note  that  T  and  M  are  mutually  exclusive. 


4.  Editing.  Up  to  three  points  may  be  replaced  in  or  added  to  the 
selected  profiles. 
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5.  Output  optims. 

Code  Meaning 

P  Plot  the  selected  profiles.  Only  one  graph 

is  generated  per  input  card;  all  the  pro¬ 
files  selected  are  plotted  on  the  same 
graph. 

N  Same  as  P,  except  that  every  3th  point  on 

the  profiles  is  annotated  with  its  profile 
index  number. 

C  Punch  the  selected  profiles. 

Each  input  card  is  a  self-contained  entity  and  specifies  the  RSVP  request 
number,  the  index  of  (me  or  more  profiles  to  be  processed,  and  for  each  index 
specified,  either  zero,  one  or  two  of  the  processing  options.  After  all  the  indices 
and  their  corresponding  processing  opticms  appears,  the  output  options  are 
specified.  This  is  followed  by  parameters  used  in  the  editing  and  filtering,  and  a 
title  used  to  identify  the  plotter /punch  output. 

The  format  of  the  card  is  given  below.  Sample  card  input  appears  on 
succeeding  pages. 


Format 

Columns 

Variable 

Name 

Description 

12 

1-2** 

3SET 

The  RSVP  request  number  of  the  set  of  pro¬ 
files  to  be  processed.  Input  cards  must  be  in 
increasing  order  by  this  number. 

8A1 

3-10 

(LSELd), 

1=1,8) 

Processing  parameters  specified  as  an  index 
(1-9)  for  each  profile  to  be  processed,  followed 
by  the  processing  options  (T,  M,  F)  for  that 
profile,  terminating  with  the  output  options 
(P,  C)  for  the  set  of  profiles.  See  the  explana¬ 
tion  above  lor  more  details.  See  the  following 
pages  for  sample  input. 

12 

F7.0 

12 

F7.0 

12 

F7.0 

11-12,** 

14-20,* 

21-22,** 

24-30,* 

31-32,** 

34-40* 

(MD(I), 

5M(1), 

1=1,3) 

-  -1 

Data  to  be  used  in  editing  sound  speed  pro¬ 
files.  Expressed  as  pairs  of  sequence  number 
within  profile  and  replacement  sound  speed  in 
meters/sec.  The  specified  editing  is  done  to 
all  profiles  selected  in  columns  3-10  above.  A 
maximum  of  3  points  can  be  replaced  or 
added.  MO(l)  cannot  be  larger  than  35.  If  no 
editing  is  desired,  leave  these  fields  blank. 

*The  value  may  appear  anywhere  in  the  indicated  columns  with  the  decimal  point 
punched. 


**The  value  must  be  right-justified  in  the  indicated  columns  with  no  decimal 
point. 


Format 

Columns 

Variable 

Name 

Description 

12 

41-42** 

LC 

Number  of  points  to  be  selected  from  each  of 
the  specified  profiles  when  the  selected 
profile  is  filtered.  Maximum  is  35.  Default  is 

15. 

F3.1 

43-45* 

TL 

Tolerance  in  meters  per  second  used  for  filter¬ 
ing  option  (F).  Maximum  is  3.75.  Default  is 

1.5.  See  description  of  Subroutine  FRITZ  in 
Section  IV.5  for  a  detailed  explanation  of  this 
parameter. 

3A10 

46-75 

NOTE2 

Up  to  30  alphanumeric  characters  used  as 
annotation  on  the  plot  and  as  an  identifier  for 
the  punched  output. 

11 

76 

NEXT 

Index  number  of  the  deep  level  to  which 
profiles  are  to  be  extended,  using  archival 
deep  sound  speeds  from  the  Type  4  record 
described  in  Section  U.l.  Next  is  used  when 
processing  option  T  is  chosen  in  columns  3-10. 
Deep  level  1  is  1000  meters,  level  2  is  2000 
meters,  and  so  on  to  level  9  at  9000  meters.  If 
NEXT  is  blank  or  0,  then  the  deep  level  is  set 
to  10000  meters. 

13 

- ,  i 

78-80 

DM  AX 

Maximum  depth,  in  hundreds  of  meters,  for 
plotting  sound  speed  profiles.  If  DMAX  is 
blank  or  0,  then  a  default  of  6000  meters  is 
assumed.  If  DMAX  is  set  to  a  negative 
number,  then  the  maximum  depth  for  plotting 
is  computed  by  the  program. 

This  parameter  card  is  repeated  for  each  set  of  data  to  be  selected,  up  to  a 
maximum  of  16  cards.  Cards  must  be  in  increasing  order  by  the  request  number 
in  columns  1  and  2. 


*The  value  may  appear  anywhere  in  the  indicated  columns  with  the  decimal  point 
punched. 

•«The  value  must  be  right-justified  in  the  indicated  columns  with  no  decimal 
point. 
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ni.  OUTPUT 

111.1  Printer  Output 

111.1.1  Standard  Output 

A  listing  of  each  profile  selected  is  printed  along  with  sequence  numbers 
for  eadi  point  in  the  profile.  A  heading  at  the  top  of  eac^t  profile  indicates  the 
original  RSVP  request  number  and  the  processing  done  on  the  profile.  Points  on 
the  profiles  which  have  been  filtered  out  are  printed  as  zero. 

See  the  following  page  for  a  sample  printout. 


III.  1.2  Error  Messages 

An  error  in  entering  the  RSVP  request  number  on  an  input  card  to 
RSVPPLP  produces  the  following  message: 

***#»*#***bad  request  no.  j 

input  card 

ABOVE  INPUT  CARD  RE3ECTED 

where  j  is  the  request  number  from  the  input  card.  The  program  continues  with 
the  next  input  card. 

Printed  by:  Main  program 

Solution:  Check  columns  1-2  to  be  sure  that  the  request  number  is  between  1 

and  16  inclusive. 


An  error  in  the  sequence  of  the  input  cards  for  RSVPPLP  produces  the 
following  message: 

♦********#rEQUEST  no.  j  OUT  OF  ORDER 
input  card 

ABOVE  INPUT  CARD  REJECTED 

where  j  is  the  request  number  from  the  input  card.  The  program  continues  with 
the  next  input  card. 


Printed  by:  Main  program 

Solution:  Check  columns  1  and  2  of  each  input  to  make  sure  that  they  are  in 

order  by  request  number. 
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FIGURE  ni-l:  SAMPLE  PRINTER  OUTPUT 


An  error  in  entering  the  options  on  an  input  card  to  RSVPPLP  produces  the 
following  message: 

♦****»»**#bad  option  o 

input  card 

ABOVE  INPUT  CARD  REJECTED 

where  o  is  the  option  from  the  input  card.  The  program  continues  with  the  next 
input  card. 

Printed  by:  Main  program 

Solution:  Check  columns  3-10  to  make  sure  that  the  options  are  M,  T,  F,  C, 

P,  or  N. 


An  error  in  the  RSVP  data  file  produces  the  following  message: 

******»«**BAD  RSVP  DATA 
RSVP  record 

The  program  aborts  after  printing  the  message  and  the  bad  record. 

Printed  by:  Main  program 

Solution:  Check  that  the  preceding  rvm  of  RSVP  was  executed  properly,  and 

that  the  correct  file  name  was  used  for  the  RSVPPLP  run. 


When  the  end  of  data  is  reached  on  the  RSVP  data  file  before  all  the  input 
cards  to  RSVPPLP  have  been  processed,  the  following  message  is  printed: 

**»ft»«****NO  MORE  DATA  ON  RSVP  RLE 

ALL  INPUT  CARDS  WERE  NOT  PROCESSED 

The  program  terminates  normally. 

Printed  by:  Main  program 

Solution:  Check  the  input  cards  for  a  request  number  that  does  not  exist  in 

the  preceding  RSVP  run. 


If  any  of  the  sound  speeds  to  be  plotted  are  outside  the  range  for  plotting, 
the  following  error  message  is  printed: 

******»#**SOUND  SPEED  POINT  i,  d,  s, 

FOR  PROFILE  INDEX  j  OUTSIDE  PLOT  BOUNDS 

where  i  is  the  sequence  number  of  the  point  within  the  profile, 
d  is  the  depth  of  the  point  in  meters, 
s  is  the  sound  speed  of  the  point  in  meters/second, 
j  is  the  profile  index  number. 

Printed  by:  Subroutine  PLPROF 

Solution:  Check  output  of  RSVP  program  for  points  with  speeds  less  than 

1200  or  greater  than  ISOO  meters/sec. 


I 


Plotted  output  is  optional;  it  occurs  only  if  the  code  P  or  the  code  N  is 
punched  in  the  appropriate  place  in  columns  3>i0  of  the  input  parameter  card. 
Only  one  graph  can  be  generated  by  each  input  card.  All  the  profiles  selected  on 
the  input  card  are  plotted  on  the  same  graph. 


See  the  descriptions  of  the  plotting  std>routines  PLOT1T,  ANNOT,  AXES, 
LEGEND,  and  PLPROF  for  details  about  the  plotted  output. 


See  the  following  page  for  a  sample  of  plotted  output. 
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ni.3  Card  Output 

Punched  card  output  is  optional.  Card  images  are  written  on  logical 
unit  3  if  the  code  C  is  punched  in  the  appropriate  place  in  columns  3-10  of  the 
input  parameter  card. 

For  each  profile,  a  set  of  card  images  is  output  consisting  of  a  header  card 
(Type  1),  a  card  containing  the  number  of  points  in  the  profile  (Type  2)  followed 
by  cards  containing  the  profile  expressed  as  depth,  sound  speed  pairs,  four  pairs 
per  card  (Type  3). 

The  format  of  the  output  cards  is  described  below.  A  listing  of  a 
sample  punched  card  output  is  shown  on  a  succeeding  page. 


Card 

Type 

Format 

Columns 

Variable 

Name 

Description 

1 

4A10 

1-40 

(ID(I,ISET), 

1=1,4) 

First  40  characters  of  the  iden¬ 
tifying  title  from  the  original 
RSVP  request  sequence  number 
ISET. 

3A10 

41-70 

(NOTE2(I), 

1=1,3) 

Alphanumeric  identification 

from  columns  46-75  of  the 
input  parameter  card. 

AlO 

71-80 

LHS(K) 

Alphanumeric  description  of 
the  type  (e.g.,  mean,  model, 
etc.)  for  the  Kth  profile 
selected. 

2 

15 

1-5 

L** 

Number  of  points  in  the  profile. 
Maximum  is  35. 

There  is  one  Type  1  card  followed  by  one  Type  2  card  for  each  profile  punched. 


*«Value  is  right-justified  in  specified  columns  with  no  decimal  point  punched. 


Card 

Type 

Format 

Columns 

Variable 

Name 

Description 

3 

8F10.3 

1-10,* 

11-20,* 

• 

• 

(DT(l),STa), 

1=1, L) 

The  sound  speed  profile,  ex¬ 
pressed  as  pairs  of  (depth, 
sound  speed)  in  (meters,  me¬ 
ters/sec. 

61-70,* 

71-80* 

Cards  of  Type  3  are  repeated,  with  up  to  4  profile  points  per  card,  until  all  points 
in  the  profile  have  been  specified,  up  to  a  limit  of  33  points. 


*The  value  is  positioned  anywhere  in  the  columns  specified  with  the  decimal 
point  punched. 


MJnCL 


riiT 


I* 

O.O)  I 
7!».0,3  1 

suo.o'a  1 

7000. OC  1 

1 

0.00 

5O.00 

l<^0.00 

*00.00 

AOO.OO 

1200.00 

17^0.00 

11 

2fl 

0.00 

*10.00 

1*10.00 

400.00 

000.00 

12U0.00 

17SO.QO 

01 

O.OQ 

so. 00 

ISO. 00 
*00.00 
AOO.Ot 
1200.00 
17S0.00 

[02 

2S 

0.00 

so.  00 

I'^O.OO 

*00.00 

•toO.OO 

12u0.00 

17«.0.03 

102 

2» 

0.00 
SO. 00 
ISO. 00 
*00.00 
HVO.OO 
1200.00 
I7S0.00 


2S 

0.00 

so.  00 

ISO. 00 

*oo.oo 

MOO. 00 
1200.00 
17S0.00 


ISV.SO 
ISIS, 60 
i5n.*o 
1*2S.00 


1S30.S3 
1S16.80 
1512. TO 
ISIS.AO 
1521.20 
1527,30 
1536,60 


1538.20 

1515.50 

1512.50 
ISl^^.TO 
1520.80 
1527.10 

1536.50 


10.00 

100.00 

MuO.OO 

ICOOO.OO 


10.00 

75,00 

200.00 

500.00 

400.00 

1300.00 

2000.00 

10.00 

75.00 

200.00 

500.00 

400.00 

1300.00 

2000.00 


1524.30  20.06  152*. 6o 

1513.30  260.00  1564.30 

1517.60  1066.06  1521.30 

1682.00  ^  ^ 

HE*»l  *•-’  6.  0. 


15*0.80 

1517.70 

1512.50 
1516.10 
1521.60 

1527.50 

1536.70 


15*2.10 

1518.60 

1513.60 
1516.20 

1521.60 

1527.60 
1536.70 


30.00  15?7.*U 
?»0.00  1506.10 
JOU0.06  IS55.60 


1537.10 
1513.*0 
1513.30 
1516.70 
1522. SO 
1525.00 
1561.00 


1536.70 

1512.50 

1511.30 

1516.50 
1520.80 

1627.30 
1536,16 


20.00 
100.00 
250.00 
600.00 
1000.00 
1*00.00 
I  2SOO.OO 
•»EaN  ♦«-  S. 


1533.60 

1512.20 
1513.00 
1516.00 

1522.20 
1528.80 
1561.00 


10.00 

75.00 

200.00 

500.00 

500.00 

1300.00 

2000.00 


10.00 

75.00 

200.00 

500.00 

400.00 

1300.00 

2000.00 


10.00 

75,00 

200.00 

500.00 

500.00 

1300.00 

2000.00 


1565.50 

1523.50 
1516.10 
1516.60 

1523.50 
1527,40 
1537.30 


20.00 

100.00 

250.00 

600.00 

1000.00 

1*00.00 

2500.00 


**£*M  6»  0. 


1526.60 

1512.70 
1513.80 
1518.10 
152*. 10 

1530.70 
15<*5.50 

D. 

1526.20 

1512.20 
1513.30 
1517.60 
1523.80 
1530.50 
15*5.50 


30.00 

125.00 

300.00 

7U0.00 

1100.00 

ISUO.OO 

JOUO.OO 


15*0.60 

1516.60 

1513.60 
1517. *0 
1522.80 
1529,20 
1561.00 


20.00 

100.00 

250.00 

6V0.0i> 

1000.00 

1*00.00 

2500.00 


M0  7ei.«MJN.M*X 


1533.00 

1513.20 

1516.30 

1518.60 

1526.60 
1530.90 
15*5.50 


15*0.60 
151*. 00 
1513.70 
1517.10 
1522.80 
1529,20 
15*1.00 


20.00 

100.00 

25U.00 

600.00 

1000.00 

J*0D.0D 

2500.00 


•*Or>El.»HlN»M*X 


1521.50 

1511.30 
1511 #90 
1516,20 

1522.30 
1528.40 

1560.30 


20.00 

100.00 

250.00 

600.00 

1000.00 

1*00.00 

2500.00 


1521.90 
Ul?.60 
1516. 3U 
1514.60 
1525.70 
1532.30 
I  1558. OU 
MEAN-SO 


30.00 

125.00 

300.00 

700.00 

1)00.00 

1500.00 

iOOO.OO 


1536.30 

1513.60 

1516.50 

1518.50 

1526.60 
1530.80 
15*9.50 


MOOEC**3I»1*8AA 


10.00 

75.00 

200.00 

500.00 

400.06 

1300.80 

2000.00 


1565,30 

1519,90 

1519.20 

1520.70 

1526.60 

1529.60 


20.00 

100.00 

250.00 

600.00 

1000.80 

1*80.80 


1519.30 

1512.30 

1513.60 

1519.10 

1525.60 

1532.10 
I  1558.00 
he AN. SO 


30.00 

12S.00 

300.00 

700.00 

1100.00 

1500.00 

3000.00 


1517.10 

1509.60 
1512.90 

1517.70 

1523.70 

1530.60 
15*8.90 


152*. SO 
1512.90 
ISIS. 00 
1520.10 
1525.80 
1532.50 
1558.00 


model 


30.00 

125.00 

300.00 

700.00 

tlOO.OO 

1500.00 

3000.00 


30.00 

125.00 

300.00 

700.00 

1100.00 

1500.00 

3000.00 


1528.30 

1513.30 
1515.20 
1519.00 
1526,00 
1532.50 
1557.60 


1516.80 

1510.10 

1513.10 
1519.30 
1525.70 
1532.00 
1557.80 


1561,30  2500.80 


15*5.90 

1518.90 

1518.70 

1521.50 

IS2S.8U 

1531.30 

15*9.90 
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3UO.OO 
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1537.60 
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153*. 10 
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FIGURE  ni-3;  SAMPLE  CARD  OUTPUT 


IV.  PROGRAM  DESCRIPTION 


IV. 1  Main  Program 

Function:  To  extend,  filter,  edit,  plot  and/or  punch  sound  speed  profiles  that 

have  been  previously  retrieved  (using  program  RSVP)  from  the 
NODC/LRAPP  Gcofile. 

COMMON  COMMON  storage  is  not  used  in  the  RSVPPLP  main  program.  It 
Storage:  is  used  only  in  subroutine  PLOTIT  lor  communication  with  sub¬ 

sidiary  plotting  routines. 

Description:  Program  RSVPPLP  reads  an  input  parameter  card.  It  checks  the  card 
for  validity,  and  then  examines  the  options  specified  in  the  card 
for  selection,  processing,  and  output  of  sound  speed  profiles. 

The  specified  profiles  are  read  from  the  data  file  output  by 
program  RSVP,  and  moved  to  a  processing  array  where  they  are 
edited  or  extended  (in  the  main  program)  or  filtered  (by  subrou¬ 
tine  FRITZ).  The  profiles,  both  before  and  after  each  stage  of 
processing,  are  stored  in  an  array  for  later  printing  using  subrou¬ 
tine  PAGE. 

U  one  of  the  output  options  on  the  input  card  is  plotting,  the 
processed  profiles  are  plotted  using  subroutine  PLOTIT  and  its 
subsidiary  plotting  subroutines  ANNOT,  AXES,  LEGEND, 
PLPROF,  and  SCALDS.  If  one  of  the  output  options  is  punched 
cards,  the  card  images  of  the  processed  profiles  are  written  to  a 
disk  file  using  subroutine  PNCH. 

Succeeding  input  cards  are  read  and  the  RSVP  profile  data 
specified  by  them  is  read  and  processed  until  there  are  no  more 
input  cards  to  be  read. 

See  the  next  page  for  a  diagram  of  the  main  program  flow. 


Initialize  counters 
and  indicators. 


Read  the  rest  of  the 
RSVP  request  data. 


RSVPPLP 
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IV.2  Subroutine  ABRT  (NAME) 

Function: 

To  abort  the  job  when  an  error  has  occurred,  after  printing  a 
message  giving  the  name  of  the  program  moidule  that  was  in 
control  when  the  error  was  detected. 

Where 

Called: 

Main  program 

Parameter 

Input: 

NAME  -  Alphanumeric  name  of  the  program  calling  ABRT. 

COMMON 

Input: 

None 

Card  Input: 

None 

File  Input: 

None 

Parameter 

Output: 

None 

COMMON 

Output: 

None 

Printer 

Output: 

The  following  message  is  printed: 

**##*#****RSVPP1,P  ABORTED  IN  name. 

File  Output: 

None 

Description: 

ABRT  prints  the  error  message  and  then  calls  a  dummy  subroutine 
ABORT  that,  on  the  CDC  6000  series  computers,  causes  the  run 
to  terminate  abr»rmally. 

Error 

Messages:  None 


IV-9 


IV.3  Subroutine  ANNOT  (ICAPl.  1CAP2.  LABEL.  ICN) 

Function:  To  put  an  annotation  at  the  bottom  of  a  graph  of  RSVP  sound 

spe^  profiles,  identifying  the  profiles  and  explaining  the  process¬ 
ing  that  has  been  done  on  each  of  the  profiles  prior  to  plotting. 

Where 

Called:  Subroutine  PLOTIT. 

Parameter  ICAPl  -  A  40-character  alphanumeric  identifier  for  the  original 

Input:  RSVP  run  that  selected  the  profiles  according  to  location 

and  month  of  the  year. 

ICAP2  -  A  30-character  alphanumeric  identifier  from  the 
RSVPPLP  card  input. 

LABEL  -  An  array  containing  up  to  5  descriptors  for  each  profile 
plotted.  Word  1  contains  the  profile  index  number  (1-9) 
used  to  annotate  each  individual  profile.  Words  2-3 
contain  the  codes  (M,T,F,E)  used  to  indicate  the  proces¬ 
sing  that  has  been  done  on  the  profile.  See  Section  IV.6 
(Subroutine  LEGEND)  for  an  explanation  of  the  profile 
index  numbers  and  the  processing  codes. 

ICN  •  The  number  of  sound  speed  profiles. 

COMMON 

Input:  None 

Card  Input:  None 

File  Input:  None 

Parameter 

Output:  None 

COMMON 

Output:  None 

Printer 

Output:  None 

Plotter  At  the  bottom  of  the  graph  of  a  set  of  profiles,  the  alphanumeric 

Output:  identifiers  in  arrays  ICAPl  and  1CAP2  are  plotted.  Under  these 

are  plotted  identifying  labels  for  each  of  the  profiles  on  the 
graph.  They  are  of  the  form: 

I  -  p  ••••  p 

1  n 

where  I  is  the  profile  index  number  used  to  annotate  the  profile, 
and  the  Ps  indicate  the  processing  performed  on  the  profile. 


IV-li 


Description:  The  subroutine  plots  the  first  caption,  followed  by  the  second 
caption,  and  then  plots  a  label  for  each  profile  on  tfie  graph. 


Error 

Messages:  None 


Draw  first  caption. 


^NQI 
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- 

1 

s 

IV.4  Subroutine  AXES 

;  Function; 

To  draw  and  annotate  depth  and  sound  speed  axes  for  RSVP  plots. 

Where 

Called: 

Subroutine  PLOTIT. 

Parameter 

Input: 

None 

COMMON 

Input: 

/DPLOT/ 

ISCALD 

-  Index  into  the  depth  scaling  arrays  used 
for  the  current  graph. 

LENY 

-  Number  of  divisions  on  the  depth  axis. 

1 

YINC 

-  Length  in  inches  of  a  division  on  the  depth 
axis. 

DPMAX 

•  Maximum  depth  in  meters  on  the  depth 
axis. 

YFAC 

-  Multiplicative  factor  to  convert  axis  tic 
annotation  to  meters. 

i 

1 

/SPLOT/ 

ISCALS 

-  Index  into  the  sound  speed  scaling  arrays 
used  for  the  current  graph. 

LENX 

*  Number  of  divisions  on  the  sound  speed 
axis. 

XINC 

>  Length  in  inches  of  a  division  on  the  sound 
speed  axis. 

SSINC 

-  Length  in  meters/sec.  represented  by  a 
division  on  the  sound  speed  axis. 

SSMIN 

-  Minimum  sound  speed  in  meters/sec.  on 
the  sound  speed  axis. 

Card  Input: 

Nwie 

File  Input: 

None 

Parameter 

Output: 

None 

COMMON 

Output: 

None 

Printer 

Output: 

None 

. . 

IV-l# 

-  J 

Plotter 

Output: 

Pile  Output: 
Description: 


Error 

Messages: 


A  depth  axis  is  drawn  with  an  origin  at  the  top  left  corner  of  the 
graph,  with  depth  increasing  downward.  This  axis  is  repeated  on 
the  right  side  of  the  graph.  Sound  speed  axes  are  drawn  on  the 
top  and  bottom  of  the  graph.  The  axes  are  labelled. 

None 

Using  the  index  ISCALD,  the  appropriate  depth  scaling  param¬ 
eters  are  picked  out  of  the  scaling  tables,  and  the  depth  axis  is 
drawn  and  numbered  on  the  left  side  of  the  graph.  Using  the 
index  ISCALS,  the  appropriate  sound  speed  scaling  parameters  are 
picked  out  of  the  scaling  tables,  and  the  sound  speed  axis  is  drawn 
and  numbered  along  the  top  of  the  graph.  Then  a  repeat  of  the 
depth  axis  is  drawn  on  the  right  si^,  and  a  repeat  of  the  sound 
speed  axis  is  drawn  on  the  bottom.  Finally,  titles  are  drawn  on 
the  left  hand  and  top  axes. 


None 
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IV.5  Subroutine  FRITZ  (D,  S,  N,  NSEL,  SR,  NR,  ER) 

Function:  To  filter  a  sound  speed  profile  by  eliminating  points  from  it  while 

still  leaving  a  good  representation  of  the  curve. 

Where 

Called:  Main  program 

Parameter 


Input: 

D 

-  An  array  containing  the  profile  depths  in  meters. 

S 

-  An  array  containing  the  profile  sound  speeds  in  me¬ 
ters/sec. 

N 

-  The  number  of  points  in  the  profile. 

NSEL 

-  The  maximum  number  of  points  that  remain  in  the 
profile  after  the  elimination  process. 

ER 

-  The  tolerance  limit  in  meters/sec.  used  in  the  elimina¬ 
tion  process. 

COMMON 

Input: 

None 

Card 

Input: 

None 

File  Input: 

None 

Parameter 

SR 

-  An  array  containing  the  sound  speeds  in  meters/sec. 

Output: 

for  the  filtered  profile.  Points  that  have  been  elimi¬ 
nated  from  the  original  array  (S)  are  replaced  by  zeros 
in  the  SR  array. 

NR 

-  The  number  of  non-zero  points  in  the  filtered  profile. 

COMMON 

Output: 

None 

Printer 

Output: 

None 

File  Output: 

None 

Description:  The  subroutine  first  checks  to  see  that  all  the  sound  speeds  are 
within  an  acceptable  range  (1410.0  -  1720.0  m/s).  Those  that  are 
outside  the  range  are  discarded.  If  there  are  less  than  3  points 
left,  the  filtering  is  not  done. 


IV-17 


The  filtering  process  goes  as  follows:  A  line  is  drawn  from  the 
first  to  the  third  point  in  the  profile,  skipping  the  second  point. 
Knowing  the  depth  of  the  second  we  can  find  where  the  second 
point  is  on  this  Une  and  thus  know  the  sound  speed  at  that  point. 
We  then  compare  this  calculated  sound  ^eed  with  the  observed 
sound  speed.  If  they  are  within  a  certain  tolerance,  then  a  line  is 
drawn  from  point  1  to  point  4.  The  second  and  tfiird  points  are 
then  checked  to  see  if  they  fall  within  the  tolerance.  If  they  do, 
then  a  line  is  drawn  from  point  1  to  point  5  and  so  on  until  a 
preceding  point  does  not  fall  within  the  tolerance. 

When  a  point  does  fall  outside  the  tolerance,  the  line  is  then 
drawn  from  the  first  point  to  the  last  point  before  the  tolerance 
was  exceeded  and  all  of  the  preceding  points  (except  the  first)  are 
eliminated.  The  endpoint  of  the  line  segment  is  then  taken  as  the 
first  point.  A  line  is  drawn  to  the  third  point  after  the  end  point 
and  the  process  is  repeated. 

Certain  points  are  left  in  the  filtered  profile,  no  matter  whether 
they  fall  in  the  tolerance  or  not.  These  points  are  the  first  and 
last  points  in  the  profile  and  the  point  with  the  minimum  sound 
speed. 


Error 

Messages:  None 
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Select  the  1-th  point. 

Add  1  to  the  number  of  selected  points. 


IV.6  Subroutine  LEGEND 


Function: 

To  put  explanatory  legends  at  the  beginning  of  a  set  of  RSVP 
profile  plots. 

Where 

Called: 

Main  program 

Parameter 

Input: 

None 

COMMON 

Input: 

None 

Card  Input: 

None 

File  Input: 

None 

Parameter 

Outputs: 

None 

COMMON 

Output: 

None 

Printer 

Output: 

None 

Plotter 

a)  A  legend  containing  an  explanation  of  the  numbers  used  to 

Output: 

annotate  each  profile  as  follows: 

1  -  MODEL 

2  -  MEAN 

3  -  MINIMUM 
-  MAXIMUM 

5  -  MEAN  -  STANDARD  DEVIATION 

6  -  MEAN  +  STANDARD  DEVIATION 

7  -  REPRESENTATIVE  1 

8  -  REPRESENTATIVE  2 

9  -  REPRESENTATIVE  3 

b)  A  legend  containing  an  explanation  of  the  codes  used  on  the 
annotation  at  the  bottom  of  the  graph  to  describe  the  processing 
done  on  each  profile: 

M  -  EXTENDED  WITH  MEAN 
T  -  EXTENDED  WITH  MEAN  AND  ARCHIVAL 
F  -  RLTERED 
E  -  EDITED 


File  Output:  None 
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Description: 

Error 

Message: 


The  two  legends  described  wider  Plotter  Output  are  drawn  one 
above  the  other. 

None 
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Page  1  of  1 


IV.7  Subroutine  PAGE  (IPG.  X,  ISET,  LBX) 


Function: 

Where 

Called: 

Parameter 

Input: 


COMMON 

Input: 

Card  Input: 

File  Input: 

Parameter 

Output: 

COMMON 

Output: 

Printer 

Output: 

File  Output: 

Description: 


To  store  an  RSVP  profile  for  printing,  and  later  to  print  it. 


Main  program 

IPG  -  An  indicator  that  is  set  1^  the  calling  program  to  0 
when  a  sound  speed  profile  is  to  be  stored  in  an 
internal  print  array,  and  to  1  when  the  sound  speed 
profiles  previously  accumulated  in  the  print  array  are 
to  be  printed. 

X  -  An  array  containing  sound  speeds  in  meters/sec.  to  be 

stored  in  the  print  array. 

ISET  •  The  RSVP  request  number  corresponding  to  the  array. 

LBX  -  An  alphanumeric  label  identifying  the  kind  of  proc¬ 

essing  that  has  been  done  on  the  profile  (e.g.  exten¬ 
sion,  filtering). 


None 

None 

None 


None 


None 

A  page  containing  up  to  12  sound  speed  profiles. 


None 

Each  time  that  subroutine  PAGE  is  called,  the  profile  in  array  X 
Is  transferred  to  an  internal  array  PRL.  An  internal  print  column 
counter,  ICOL,  is  incremented  by  1  for  each  array  stored.  When 
the  number  of  columns  reaches  12,  or  when  the  parameter  IPG  is 
set  to  1,  a  page  of  sound  speed  profiles  is  printed,  one  profiie  pw 
column,  with  the  request  s^uence  number  QSET),  and  the  identi¬ 
fying  label  (LBX)  at  the  top  of  each  column. 


Error 

Messages:  None 


Increment  print 
column  number. 


Print  a  page  consisting 
of  a11  remaining  data 
for  this  run. 


IV.8  Subroutine  PLDASH  (XI.  Yl.  X2.  Y2) 


Function:  To  plot  a  dashed  straight  line  from  one  point  to  another. 

Parameter 

Input:  XI  -  The  X-coordinate  in  inches  of  the  beginning  point. 

Y 1  -  The  Y-coordinate  in  inches  of  the  beginning  point. 

X2  -  The  X-coordinate  in  inches  of  the  ending  point. 

Y2  -  The  Y-coordinate  in  inches  of  the  ending  point. 

COMMON 
Input:  None 

Card  Input:  None 

File  Input:  None 

Parameter 
Output:  None 

COMMON 
Output:  None 

Printer 

Output:  None 

Plotter 

Output:  A  line  of  dashes  between  (XI,  Yl)  and  {X2,  Y2). 

File  Output:  None 

Description:  The  length  of  a  dash  and  the  spacing  between  dashes  au-e  constants 
in  the  program.  The  number  of  dashes  and  the  X  and  Y  distances 
between  ^shes  are  calculated  using  these  constants.  Then  the 
dashed  line  is  drawn. 
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Move  with  pen  up  to 
the  end  point  of  the 
line. 


IV.9  Subroutine  PLDOT  (XI,  Yl,  X2,  Y2) 


Function: 

To  plot  a  dotted  straight  line  from  one  point  to  another. 

Parameter 

Input: 

XI 

-  The  X-coordinate  in  inches  of  the  begirviing  point. 

Yl 

-  The  Y-coordinate  in  inches  of  the  begiruiing  point. 

X2 

-  The  X-coordinate  in  inches  of  the  ending  point. 

Y2 

-  The  Y-coordinate  in  inches  of  the  ending  point. 

COMMON 

Input: 

None 

Card  Input: 

None 

File  Input; 

None 

Parameter 

Output: 

None 

COMMON 

Output: 

None 

Printer 

Output: 

None 

Plotter 

Output;  A  line  of  dots  between  (XI,  Yl)  and  (X2,  Y2). 

File  Output;  None 

Description:  The  distance  between  dots  is  a  constant  in  the  program.  Using 
this  constant,  the  number  of  dots  and  the  X  and  Y  distances 
between  dots  are  calculated.  Then  the  dotted  lir>e  is  drawn. 
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IV.IO  Subroutine  PLOTIT  (D.  SS,  ICN,  ICAPl,  ICAP2,  LABEL,  LOPT,  DMAX) 


Function: 

Where 

Called: 

Parameter 

Input: 


Cvd  Input: 

File  Input: 

Parameter 

Output: 


To  plot  a  set  of  RSVP  sound  speed  profiles  on  one  graph. 


Main  program 


D  -  An  array  containing  the  profile  depths  in  meters. 
These  depths  are  the  same  for  all  profiles. 

SS  -  An  array  containing  sound  speeds  in  meters/sec.  for 
all  profiles.  A  zero  value  indicates  that  there  is  no 
sound  speed  data  at  that  point. 

ICN  -  The  number  of  sound  speed  profiles  stored  in  array  SS. 
Maximum  value  is  10. 

ICAPl  '  A  40-character  alphanumeric  identifier  for  the  RSVP 
run  that  selected  the  profiles  according  to  location 
and  month  of  the  year. 

ICAP2  -  A  30-character  alphanumeric  identifier  from  the 
RSVPPLP  card  input. 

LABEL  -  An  array  containing  up  to  5  descriptors  for  each 
profile  plotted.  Word  1  contains  the  profile  index 
number  (1-9)  used  to  annotate  each  individual  profile. 
Words  2-5  contain  the  codes  (M,  T,  F,  E)  used  to 
indicate  the  processing  that  has  been  done  on  the 
profile.  See  Section  IV.6  (Subroutine  LEGEND)  for  an 
explanation  of  the  profile  index  numbers  and  the 
processing  codes. 

LOPT  -  A  user -specified  plot  option  indicator.  If  LOPT  is  set 
to  P,  the  profile  is  not  annotated.  If  LOPT  is  set  to  N, 
the  profile  is  annotated  at  every  5th  point  with  the 
profile  index  number. 

DMAX  -  The  user-specified  maximum  depth  in  meters  for  plot¬ 
ting.  If  DMAX  is  set  to  0,  a  default  of  6000  meters  is 
assumed.  If  DMAX  is  set  negative,  the  maximum 
depth  is  determined  from  the  data  in  the  D  array. 

None 

None 


DMAX  -  The  maximum  depth  in  meters  for  plotting. 
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/ 


COMMON 

Output: 


Subroutine  PLOTIT  sets  up  COMMON  storage  tables  that  are  used 
by  its  subsidiary  subroutines  ANNOT,  AXES,  LEGEND,  PLPROF, 
and  SCALDS  to  choose  scaling  factors  for  the  plotting  of  depth 
and  sound  speed.  There  is  one  labelled  COMMON  block  for  depth, 
and  another  for  sound  speed. 


Depth  Scaling 

/DPLOT/<  NOX  -  Number  of  available  depth  scales.  Cur¬ 
rently  set  to  6. 

ISCALD  -  Index  into  the  depth  scaling  arrays  used 
for  the  current  graph. 

The  following  arrays  contain  NOX  entries,  one  for  each  depth 
scale: 


LENY  -  Numbn-  of  divisions  on  the  depth  axis. 

YINC  -  Length  in  inches  of  a  division  on  the  depth 

axis. 

DPINC  -  Length  in  meters  represented  by  a  divi¬ 
sion  on  the  depth  axis. 

DPMAX  -  Maximum  depth  in  meters  on  the  depth 
axis. 

YFAC  -  Multiplicative  factor  to  convert  axis  tic 
annotation  to  meters. 

The  values  in  these  arrays  are  set  as  follows: 


Array 

Index 

LENY 

YINC 

DPINC 

DPMAX 

YFAC 

1 

6 

1.4 

40. 

240. 

I. 

2 

6 

1.4 

100. 

600. 

100. 

3 

6 

1.4 

200. 

1200. 

100. 

4 

6 

1.4 

300. 

3000. 

100. 

5 

6 

1.4 

1000. 

6000. 

1000. 

6 

6 

1.4 

2000. 

12000. 

1000. 
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Sound  speed  scaling 


/SPLOT/  NOY 

•  Number  of  available  sound  speed  scales. 
Currently  set  to  6. 

ISCALS 

-  Index  into  the  sound  speed  scaling  arrays 
used  for  the  current  graph. 

The  following  arrays  contain  NOY  entries,  one  for  each  sound 
speed  scale: 

LENX 

-  Number  of  divisions  on  the  sound  speed 
axis. 

XINC 

-  Length  in  inches  of  a  division  on  the  sound 
speed  axis. 

SSINC 

-  Length  in  meters/sec.  represented  by  a 
division  cm  the  sound  speed  axis. 

S5MIN 

-  Minimum  sound  speed  in  meters/sec.  on 
the  sound  speed  axis. 

SSMAX 

-  Maximum  sound  speed  in  meters/sec.  on 
the  sound  speed  axis. 

The  values  in  these  arrays  are  set  as  follows: 


Array 

Index 

LENX 

XINC 

SSINC 

SSMIN 

SSMAX 

1 

6 

1.0 

20. 

1440. 

1560. 

2 

6 

1.0 

20. 

1480. 

1600. 

3 

6 

1.0 

40. 

1400. 

1640. 

6 

1.0 

40. 

1440. 

1680. 

3 

6 

1.0 

40. 

1480. 

1720. 

6 

6 

1.0 

100. 

1200. 

1800. 

Printer 

Output:  None 

Plotter 

Output:  Output  from  PLOTIT  is  one  graph  of  sound  speed  profiles.  The 

axes  are  annotated,  and  identifying  titles  are  drawn.  More  than 
one  sound  speed  profile  may  be  drawn  on  the  graph.  Each  profile 
may  be  annotated. 


File  Output:  None 
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Description:  Subroutine  PLOTIT  calls  subroutine  SCALDS  to  determine  the 
scaling  to  be  used  for  plotting  depth  and  sound  speed.  It  then 
calls  subroutine  AXES  to  draw  and  annotate  the  axes  for  the 
graph,  subroutine  PLPROF  to  plot  the  profiles,  and  subroutine 
ANNOT  to  draw  captions  at  the  bottom  of  the  graph. 

Error 

Messages:  None 
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IV.ll  Subroutine  PLPROF  (D.  S,  LABP,  LOPT) 


Function:  To  plot  an  RSVP  sound  speed  profile. 

Where 

Called:  Subroutine  PLOUT. 

Parameter 

Input:  D  -  An  array  containing  profile  depths  in  me¬ 

ters. 

S  -  An  array  containing  sound  speeds  in  me¬ 

ters/sec.  A  zero  value  indicates  that 
there  is  no  sound  speed  data  at  that  point. 

LABP  -  A  profile  index  number  used  for  annota¬ 
tion  on  the  profile  plot.  See  subroutine 
LEGEND  for  a  description  of  this  index 
number. 

LOPT  -  A  user -specif  led  plot  option  indicator.  If 
LOPT  is  set  to  P,  the  profile  is  not 
annotated.  If  LOPT  is  set  to  N,  the 
profile  is  annotated  with  the  profile  index 
number  (LABP)  at  every  5th  point. 

COMMON  /DPLOT/  ISCALD  -  Index  into  the  depth  scaling  arrays  used 

Input:  for  the  current  graph. 

YINC  -  Length  in  inches  of  a  division  on  the  depth 
axis. 

DPINC  -  Length  in  meters  represented  by  a  divi¬ 
sion  on  the  depth  axis. 

DPMAX  -  Maximum  depth  in  meters  on  the  depth 
axis. 

/SPLOT/  ISCALS  -  Index  into  the  sound  speed  scaling  arrays 

used  for  the  current  graph. 

XINC  >  Length  in  inches  of  a  division  on  the  sound 
spe^  axis. 

SSINC  -  Length  in  meters/sec.  represented  by  a 
division  on  the  sound  speed  axis. 

SSMIN  >  Minimum  sound  speed  in  meters/sec.  on 
the  sound  speed  axis. 

SSMAX  -  Maximum  sound  speed  in  meters/sec.  on 
the  sound  speed  axis. 
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Card  Input: 

File  Input: 

Parameter 

Output: 

COMMON 

Output: 

Printer 

Output: 

File  Output: 

Plotter 

Output: 


Description; 


Error 

Messages: 


None 

None 


None 


None 

Error  Message.  See  below  for  details. 


None 


Subroutine  PLPROF  plots  one  profile.  If  LOPT  is  set  to  N,  a  plot 
annotated  at  every  5th  point  with  the  profile  index  number 
(LABP)  is  produced.  If  LOPT  is  set  to  P,  a  plot  is  produced  with 
no  annotation. 

The  profile  points  that  have  zero  sound  speed  or  depths  deeper 
than  the  maximum  depth  for  plotting  are  eliminated  from  the 
profile.  Then,  using  the  scaling  indicated  by  the  values  in  ISCALD 
and  ISCALS,  the  profile  is  plotted,  either  annotated  or  not, 
depending  on  the  setting  of  LOPT. 

The  subroutine  checks  that  the  sound  speeds  to  be  plotted  are 
within  the  specified  minimum  and  maximum  sound  speed  values. 
If  any  are  not,  an  error  message  is  printed,  and  processing 
continues.  See  Section  III.  1.2  for  a  detailed  description  of  the 
error  messages. 


Pick  up  depth  and 
sound  speed  scaling 
from  tables. 


With  pen  down,  plot 
next  point  in  plotting 
array. 


With  pen  down,  plot 
remaining  points  in 
plotting  array. 


1V.12  Subroutine  PNCH  (D,  S,  ICAPl,  ICAP2,  ICAP3) 


Function: 


Where 

Caiied: 

Parameter 

Input: 


COMMON 

Input: 

Card  Input: 

File  Input: 

Parameter 

Output: 

COMMON 

Output: 

Printer 

Output: 

File 

Output: 

Description: 


Error 

Messages: 


To  write  a  sound  speed  profile  to  a  disk  file  in  card  image  form 
for  punching. 


Main  program. 

D  -  Array  containing  profile  depths  in  meters. 

S  -  Array  containing  profile  sound  speeds  in  meters/sec. 

ICAPl  -  A  40-character  alphanumeric  identifier,  from  the  orig¬ 
inal  RSVP  request.  Written  in  columns  1-40  of  the 
first  card  image  for  the  profile. 

ICAP2  -  A  30-character  alphanumeric  identifier,  from  the  in¬ 
put  parameter  card.  Written  in  columns  41-70  of  the 
hrst  card  image  for  the  profile. 

ICAP3  -  10  alphanumeric  characters  identifying  the  profile 
type  and  the  processing  done  on  it.  Written  in  columns 
71-80  of  the  first  card  image  for  the  profile. 


None 

None 

None 


None 


None 


None 

On  logical  unit  3,  card  images  of  the  sound  speed  profile.  See 
Section  II1-3  for  a  description  of  the  file. 

A  header  card  containing  identifiers  for  the  profile  is  written, 
followed  by  a  card  containing  the  number  of  points  in  the  profile. 
Then  the  profiles  are  written  out,  4  points  per  card  image. 


None 
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ENTER 


PNCH 
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IV.13  Subroutine  SCALDS  (D,  SS,  ICN,  DMAX) 


Function:  To  determine  the  scaling  for  depths  and  sound  speeds  on  a  multi¬ 

profile  plot. 

Where 

Called:  Subroutine  PLOTIT 

Parameter  D  -  An  array  containing  profile  depths  in  meters.  These 
IniHit:  depths  are  the  same  for  all  profiles. 

SS  -  An  array  containing  sound  speeds  in  meters/sec.  for 
all  profiles.  A  zero  value  indicates  that  there  is  no 
sound  speed  data  at  that  point. 

ICN  -  The  number  of  sound  speed  profiles  stored  in  array  SS. 
Maxirvium  value  is  10. 

DMAX  -  The  user -specified  maximum  depth  in  meters  for  plot¬ 
ting.  If  DMAX  is  set  to  zero,  a  default  of  6000  meters 
is  assumed.  If  DMAX  is  set  negative,  the  maximum 
depth  is  determined  from  the  data  in  the  D  array. 

COMMON  /DPLOT/  NOX  -  Number  of  available  depth  scales. 

Input: 

DPMAX  -  Array  of  maximum  depths  in  meters,  one 
for  each  depth  scale. 

/SPLOT/  NOY  -  Number  of  available  sound  speed  scales. 

SSMIN  -  Array  of  minimum  sound  speeds  in  me¬ 
ters/sec.,  one  for  each  sound  speed  scale. 

SSMAX  -  Array  of  maximum  sound  speeds  in  me¬ 
ters/sec.,  one  lor  each  sound  speed  scale. 

Card  Input:  None 

File  Input:  None 

Parameter 

Output:  DMAX  -  Maximum  depth  in  meters  for  plotting. 

COMMON  /DPLOT/  ISCALD  -  An  index  into  the  tables  used  for  scaling 

Output:  the  depth  for  plotting. 

/SPLOT/  ISCALS  -  An  index  into  the  tables  used  for  scaling 

the  sound  speed  for  plotting. 

Printer 

Output'  None 

File  Output:  None 
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A0-A109  855  OCEAN  DATA  SYSTEMS  INC  ROCKVILLE  MO 
RSVP  UPGRADE  AND  MAINTENANCE. (U) 
dec  81  J  WEBSTER 


UNCLASSIFIED 


N00014-7B-C-0818 

ML 


Description: 


Error 

Messages: 


The  subroutine  searchs  through  the  multi-profile  array  to  find 
minimum  and  maximum  sound  speeds  over  all  the  profiles.  It  then 
searches  through  the  tables  of  sound  speed  minima  and  maxima 
until  it  finds  a  set  between  which  these  computed  minimum  and 
maximum  will  fit.  If  a  fit  cannot  be  found,  the  broadest  scale 
available  is  selected. 

The  subroutine  finds  the  deepest  depth  at  which  sound  speed  data 
exist  in  the  multi-profile  array.  If  DMAX  has  been  specified  as  a 
negative  value,  then  DMAX  is  reset  to  this  maximum  value.  If 
DMAX  has  been  specified  as  zero,  then  DMAX  is  reset  to  6000. 
Otherwise,  DMAX  is  left  as  originally  specified.  The  value  of 
DMAX  is  then  used  to  search  through  the  tables  of  depth  maxima 
to  find  a  value  which  will  fit  the  value  of  DMAX.  11  none  is 
found,  then  the  deepest  maximum  available  is  used. 


None 
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IT 


Is  maximum  plot 
depth  ■  07 


Set  maximum  plot 
depth  to  actual 
maximum  of  the  data. 


Set  ISCALO  as  the 
pointer  into  the  table 


CROSS  REFERENCE  CHARTS 


V.2  COMMON  Cross  Reference 


R  —  The  indicated  program  refers  to  the  specified  COMMON  variable. 

D  —  The  specified  COMMON  variable  is  defined  in  the  indicated  program. 


VI.  USER’S  GUIDE 
VI.1  File  Identification 

The  FORTRAN  source  and  the  binary  code  for  RSVPPLP  are  stored  at 
DTNSRDC  and  at  Eglin.  The  file  names  are  the  same  at  both  places.  The 
file  identifiers  is  shown  below. 


FORTRAN  Source 
(Program  Library) 

Binary  Code 

ID  at 
Eglin 

SNat 

Eglin 

ID  at 
NSRDC 

PVEHRSVPPLOTPUNCHPL 

PVEHRSVPPLOTPUNCHLGO 

KM 

USETl 

PVEH 

VI.2  Sample  3obs 

Sample  |ob  streams  ior  executing  RSVPPLP  at  Eglin  and  DTNSRDC  are 
presented  on  the  following  pages. 

1.  Run  at  Egiin  to  produce  a  Zeta  plot  and  to  punch  sound  speed  profiies. 


XieT*T37S*M*  lUflHIMOTO 

Cohhcmt. 
cohhfmt. 

COMMENT. 

comment. 

COMMENT. 

AT  TSCMtLOO • P VeHRSVPPLOTPUNCHUGO  * 10>KM«SNaU$CT 1 •MH*1 . 
ATT ACM • TAPCl • PS VPOUT • 1 OaKM • SN>USCT 1 «MRs 1 . 

A  T  T ACM  t Plot  f ICT APlOT • l 0>MZ  «  SN>USeT2  «  MRa 1 . 
RC0UESTiPLTC0Sf«0. 

RC0UeST«TAPC3*«0. 


laaa  eziuvriaNzc 


RUN  RSVPPLP  AT  LGLlN  TO 
make  a  /e>a  plot  and  10 
punch  sound  speed  PMOriLES. 


L0A0*PL0T* 

LOO. 

R0UTe«TAPE3*0CaPU«TtU«A0* 

ROUTE tPLTCOSiDCaPUtTlOaAOtECaeOCOL*  . 

0000000000000000000000^ - 7/#/f 

OAitrcP 

OATTPCP 

OAATrcP 

•AOTrcP  ^  ^  .  . 

iOOOOOOOO«0«00000000004r>**4/7/«/f  €UMVt 


MODEL 

REPl 

REM? 

REPS 


tA 


2.  Run  at  Eglin  to  punch  sound  speed  profiles. 


IIOTtTjrStPA*  MASHIMOTO 

comment, 
comment, 
comment, 
comment, 

ATTACH«LOO*PVCHR5VPPLOTPUNCHLOO • lOaKHt tNaUSCTl •NRal , 
AT  T ACM  tTAPC 1 *RSyPOUT  t lOaKH  t  SNaUSCT I •MRa 1 , 
ReeuCST«TAPC)**0. 

LOO. 

ROUTE tTAPCJ»0CaPU»TIOaAO.  ^  ^  - 

••00000000000000000000^ —  7/r/f 
•AlTfC  WCL 

•ATTPC  «P» 

•ROTfC  NCR? 

....... 

••••••••••••••••••••••4^  <  AA/f 


UOR  REIMVTIOHZC 


run  RSVPPLP  AT  ICLIN  TO  PUNCH 
SOUND  SPEtO  PKOriLES. 
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3.  Run  at  l^lin  to  produce  a  deferred  plotting  file  to  preview  on  the 
Tektronix. 


Xia7*T3TS«r4i.  NAWINOTO 

COMMrNT. 

COMMCMT. 

COMHCNT. 

COMNCMT. 

COMMITUT. 

ATTACH*La0«PVEHRSVPPL0TPUNCHL60tI0«KM«SNBUSETl*HR>l. 

ATTACMtTAPEltRSVPOUTiIOBXHtSNBUSETltNRBl. 

REQUEST  *  PLTCOS  f 'PP « SN>USET 1 . 

ATTACM>TEK«TEKPROCPILE«IO>KV«SN>USET2«MRb1» 

REOINtMAKEtTCKiRSVPPLP. 

catalog. PLTC0S»RSVPPL0TTEST. IObKM.  . 

ooouoooooooooooooooouo4--~'7/r/f 

•AITPP  NOPEL 

OaTTPP  MEPl 

04STFP  NiP2 

OAQTPP  .  ,  ,  ^  HEP3 

•eoeoooooo«ooooooooouo<4 — </7/r/f  CAttC 


laOR  eZlRVI 16H2C 


RUN  RSVPPLP  AT  EGLIN  TO  GET  A 
OePERMEO  PLOTTING  FILE  POK  USE 
WITH  IHE  IEK1RON1X  PNEVIEWEK. 


<6 


1 


4.  Run  at  Carderock  to  produce  a  Calcomp  plot  and  to  punch  sound  speed 
profiles. 


RUUAV«Cna0900  *  rSO  tPA  *nti • 
CN ANOr , PU JA y 1 WGG T I G« 1 « . 
£0HHENT. 

comment. 

COMMENT. 

comment. 


run  ksvpplp  at  carderock 
TO  PRODUCE  calcvmp  Plots 
AND  PUNCMLO  SOUND  SPEED  PMOFIlES. 


ATTACHtLGOfPVEHRSVPPLUTPUNCMLOOtlOBPVEHtMHBl. 
ATT ACHt  tape 1 tRSVPOUT  f IObPVEN*  NWbi , 

LABEL«TAPE20«0BNVtR*RlNC«LBRSVPPL0TtVSNBCR0G37* 

R0UTE«TAPC3t0CsPU»0CFtT10BC. 

ATTACMtCALCQSA. 

LIBRART»CALC«3G. 

too. 

Gioooomooooooooooooo^— 7/^/f  cmL 

tllTPCP 
G22SGCP 
•JJ3ACP 
RATTfCR 
GAGTPCR 
•STTfCR 


MOREL 

MEAN  *.•  t.  D. 

nODfLtNINfMAX 

REPl 

REPE 

NCPl 

IlCPt 
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OCEAN  DATA  SYSTEMS,  INC. 

6000  EXECUTIVE  BLVD..  ROCKVILLE.  MARYLAND  208S2  •  301/881-3031 


MEMORANDUM 


To; 

From: 

Date: 


George  Kerr,  NORDA  Code  321 
Jacqueline  Webster  ^  (AT” 
December  9,  19S1 


The  Cyber  Control  Language  (CCL)  procedure  for  running  RSVP  has  been 
modified  as  you  requested.  It  is  stored  as  cycle  4  of  file  RSVPPROC,  ID=KT.  A 
second  procedure  called  GETDATA  was  adtted  to  the  procedure  file.  The  proce¬ 
dure  for  running  RSVP  (PRSVP)  invokes  the  GETDATA  procedure. 

You  no  longer  need  to  specify  the  magnetic  tape  Volume  Serial  Number 
(VSN)  on  your  BEGIN  command  when  running  the  procedure.  Instead,  you  can 
specify  the  ocean  basin  on  the  BEGIN  command,  and  the  procedure  will  automat¬ 
ically  mount  the  appropriate  tape.  If  there  is  a  problem  with  the  tape,  the 
procedure  will  mount  the  backup  tape.  You  can,  however,  still  specify  your  own 
tape  serial  number  on  the  BEGIN  card,  and  this  will  override  the  automatic 
mount. 

The  order  of  the  parameters  on  the  BEGIN  command  has  been  changed 
slightly.  The  following  is  a  description  of  the  parameters  for  the  procedure: 

.PRCX:,  PRSVP,  OCNBASN,  STORFIL,  TAPENUM, 

EXTRACT=NO,  SAVE=NO,  PLOT=NO,  CARDS=NO,  PROCnL=#FILE 

OCNBASN  -  Name  of  the  ocean  basin  -  ATLARC,  CARIB,  CENATL, 
CENPAC,  INDIAN,  MEDSEA,  NORATL,  NORPAC,  PH10K3P, 
or  SATLPAC 

STORFIL  —  Name  of  the  permanent  file  on  which  the  extracted  RSVP 
data  is  to  be  stored.  The  ID  will  be  KM.  This  parameter  is 
used  only  when  the  parameter  SAVE  is  set  to  YES. 

TAPENUM  —  Volume  Serial  Number  (VSN)  of  tape  to  be  used  as  input  to 
RSVP.  This  parameter  is  optional;  when  used,  it  overrides 
the  automatic  tape  mount. 

EXTRACT  —  Set  to  YES  to  extract  RSVP  data,  set  to  NO  otherwise. 

SAVE  —Set  to  YES  to  save  the  extracted  RSVP  data,  set  to  NO 
otherwise. 


1 


PLOT 


1 

If 


—  Set  to  YES  to  plot  the  extracted  RSVP  data,  set  to  NO  other¬ 
wise. 


CARDS  —Set  to  YES  to  select  sound  speed  profiles  from  the  extracted 
RSVP  profiles,  and  to  store  these  selected  profiles  as  card 
images  on  disk  in  a  format  accepted  by  the  Model  Interface 
program. 

Attached  are  a  tabulation  of  the  tape  serial  numbers  and  a  listing  of  the 
CCL  procedure. 
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RSVP  Input  Tapes 


Original 

Backup 

Ocean 

Tape 

Eglin 

NORDA 

Eglin 

NORDA 

Basin 

Label 

VSN 

VSN 

VSN 

VSN 

ATLARC 

PUJAATLARC 

5388 

1346 

11697 

1307 

CARIB 

PU3ACARIB 

4050 

1313 

12100 

1222 

CENATL 

PU3ACENATL 

0774 

1499 

0863 

1495 

CENPAC 

PU3ACENPAC 

0462 

1326 

11020 

1212 

INDIAN 

PU3AINDIAN 

3799 

1228 

5150 

1496 

MEDSEA 

PU3AMEDSEA 

4036 

1309 

16848 

1357 

NORATL 

PU3ANORATL 

0400 

1493 

6135 

1287 

NORPAC 

PU3ANORPAC 

10575 

1234 

16386 

1355 

PHIOK3P 

PU3APHIOK3P 

10996 

1213 

12645 

1501 

SATLPAC 

PU3ASATLPAC 

15056 

1368 

16382 

1334 
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